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Abstract 

The radiative seesaw model with an inert doublet has been shown to be attractive from 
a viewpoint of both neutrino masses and cold dark matter. However, if we apply this model 
to the explanation of the positron excess in the cosmic ray observed by PAMELA, a huge 
boost factor is required although it can be automatically explained that no anti-proton 
excess has been observed there. We consider an extension of the model to enhance the 
thermally averaged annihilation cross section without changing the features of the model 
favored by both the neutrino oscillation and the relic abundance of dark matter. It is 
shown that the data of PAMELA and Fermi-LAT can be well explained in this extended 
model. Constraints from gamma ray observations are also discussed. 
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1 Introduction 



The existence of dark matter (DM) now gives us a clear motivation to investigate physics 
beyond the standard model (SM). Although we know its relic abundance in the present 
universe PP, its nature is not clarified except that DM should be cold. However, recent 
observational data on cosmic rays may give us interesting information on its mass and 
interactions. PAMELA has reported positron excess in the cosmic ray at the 6 - 100 GeV 
range in comparison with the expected background [2]. However, it has observed no 
anti-proton excess. The preliminary report of Fermi-LAT also suggests that the total 
flux of positrons and electrons is larger than the expected background at regions of 60 - 
1000 GeV [3], although any bump shown in the ATIC result [4 J is not found in that flux. 
If we consider these results are caused by the decay or the annihilation of DM but not 
by astrophysical origins, they are expected to give us crucial information on the nature 
of DM. However, it has been pointed out that there is a difficult problem if we try to 
understand these results on the basis of DM physics. 

In case of the DM decay, DM life time should be extremely long such as O(10 26 ) sec 
in order to explain the PAMELA positron excess [5j E]. It is not so easy to answer 
how such a long lifetime can be naturally realized in the ordinary models for elementary 
particles, although there are several proposals to solve this problem. In case of the DM 
annihilation, its relic abundance requires the thermally averaged annihilation cross section 
(av) to be O(10 -26 ) cm 3 /sec at its freeze-out period where the typical DM velocity t>DM 
satisfies ^dm/c ~ 0.2. On the other hand, the positron excess found in the PAMELA 
experiment requires (av) to be O(10 -23 ) cm 3 /sec for the DM in our Galaxy where Vdm/c ~ 
10~ 3 is expected for the averaged DM velocity. This means that there should be some 
large enhancement introduced as a boost factor usually, which may be caused by particle 
physics, or astrophysics, or both of them. There have been several analyses on this point 
[7J [8] . In this paper, we focus our attention on the particle physics side and study a way 
to overcome this difficulty in a model for both neutrino mass and DM. 

The radiative seesaw model proposed in [9] gives an attractive scenario for the neutrino 
mass and mixing. They are explained by new physics at TeV scales in this model. The 
model includes DM candidates whose stability is guaranteed by a discrete symmetry. 
It forbids bare Dirac mass terms of neutrinos and then is related to the smallness of 
neutrino masses. The model has been studied from various points of view [101 HU [T2| [13]. 
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However, if we apply this model to explain the PAMELA data by the DM annihilation, 
we confront the boost factor problem, unfortunately. In this model the annihilation cross 
section has a dominant contribution from its p-wave part. Since the p-wave contribution 
becomes smaller for smaller DM relative velocity v, the situation is much worse than 
the s-wave case. In fact, this requires a huge boost factor of O(10 6 ) to explain the 
PAMELA data on the basis of this model unless there are some additional astrophysical 
effects [T2j [T3] . On the other hand, we should also remind the reader that the model 
has an interesting feature favored by the PAMELA data, that is, the DM can annihilate 
only to leptons. Moreover, if we impose constraints on the model from the lepton flavor 
violating processes like /x — > ej, should not be yielded as the primary final states of 
the annihilation. Positrons originated from this DM annihilation are generated through 
the decay of /i + and [13]. Model independent analyses suggest that this feature is 
again favored by both data of PAMELA and Fermi-LAT [8j. Thus, it is an interesting 
subject for this model to find some solutions for this boost factor problem by extending 
the model without disturbing these nice features of the modelO 

In this paper we propose a simple extension of the model, which makes the Breit- 
Wigner enhancement of thermally averaged annihilation cross section possible. In that 
model we show that both data of PAMELA and Fermi-LAT can be well explained without 
assuming an unknown huge boost factor as long as the mass of a scalar field is finely 
tuned. The enhanced annihilation cross section may also predict the large flux of gamma 
ray which is associated with the DM annihilation. We discuss the consistency with the 
data for the diffuse gamma ray from observations in the EGRET, HESS and Fermi-LAT 
experiments. 

The following parts of the paper are organized as follows. In section 2 we fix the 

radiative seesaw model and discuss how all the neutrino oscillation data, lepton flavor 

violating processes and the DM relic abundance can be consistently explained. After 

that we extend the model to enhance the DM annihilation cross section in the present 

Galaxy. In section 3 we address the features required for the explanation of the data of 

PAMELA and Fermi-LAT. We also discuss the consistency between the diffuse gamma 

ray flux expected in the model and the present experimental data. We summarize the 

d One solution has been proposed by considering the decaying DM in a supersymmetric extension of 
the model [6]. 
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paper in section 4. 

2 Breit-Wigner enhancement 
2.1 A radiative seesaw model 

The radiative seesaw model considered here is an extension of the SM with an inert 
doublet i] (an additional SU(2)i doublet scalar with (rj) = 0) and three gauge singlet 
right-handed fermions iV& (A; = 1, 2, 3) [9]. In order to forbid tree-level Dirac masses for 
neutrinos, we impose Z 2 symmetry on the model. An odd charge of this Z 2 symmetry is 
assigned to all of these new fields, although an even charge is assigned to all of the SM 
contents. Both interaction Lagrangian Cn relevant to N k and scalar potential V invariant 
under the imposed symmetry are written as 

C N = - (M k W k P R N k + M k W k P L N c k ) - (h ak Z V P R N k + h.c), 

V = mj0V + ^^ + Ai(0V) 2 + A2(^) 2 + A 3 (0V)(^) + A4(0 t ^)(^V) 

+ y[W + h.c.]. W 

where £ a and (ft stand for a left-handed lepton doublet and an ordinary Higgs doublet 
scalar, respectively. Coupling constants and masses of singlet fermions are assumed to be 
real, for simplicity. C N is assumed to be written by using the basis in which a charged 
lepton mass matrix is diagonal. 

The model has several interesting features [TUl [111 H2j EE]. First, neutrino masses are 
generated through one-loop diagrams as 

■Ma/3 = ^kh ak h/3 k , 
k 

Ak = ^Ml / (Ml) I(x) = -E-(i + (2) 

where M 2 = m 2 + (A3 + A4)(0) 2 . This neutrino mass matrix can explain the neutrino 
oscillation data well as long as we set appropriate values for the parameters A5, h ak , M k 
and M v . We note that A 5 should be very small to generate desired neutrino masses. 
However, this tuning is not so bad nature since it can be controlled by a global symmetry 
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which appears if we make A5 zerojfl Second, the lightest one of N k can be DM since 
its stability is guaranteed by the Z 2 symmetry. Its relic density can be adjusted to 
the one required by WMAP for the same parameters used for the explanation of the 
neutrino oscillation data. Third, these are also consistent with the constraints imposed 
by the lepton flavor violating processes such as \i — > and r — > /ry, if neutrino Yukawa 
couplings h ak have certain flavor structure^ 

In order to show these aspects concretely, we consider an example of such flavor 
structure for neutrino Yukawa couplings as 



hei = 0, hai = h Ti (i = 1,2); h 



e3 



-h 



t-3) 



h 



//3 



-h 



t-3- 



(3) 



In this case the neutrino mass matrix can be written as 

/ \ ( I \ 



M 



1 1 

Vo 1 1 j 



(h 2 Tl At + h 2 r2 A 2 ) + 



1 



1 1 -1 

V-i -1 1 / 



h 2 A 



3 JV 3, 



(4) 



and the tri-bimaximal neutrino mixing is automatically realized for the neutrino mass 
matrix (T5]) [13]. Moreover, only two mass eigenvalues take nonzero values. Thus, the 
neutrino oscillation data can be consistently explained as long as the following conditions 
are satisfied: 



h 2 Tl A x + h\ 2 A 2 ~ 2.5 x 1(T 2 eV, h 2 r3 A 3 ~ 2.9 x 10~ 3 eV. 



(5) 



These come from the required values for Am^ tm and Amg 0lar , respectively. We need 
to consider the constraints from both the lepton flavor violating processes and the DM 
relic abundance under these conditions. The relation of Yukawa couplings h Q k to other 
parameters A5, and M v is also determined through these constraints. When we apply 
eq. (Q to the analysis, it may be useful to note that these give the constraints on the 
value of h 2 . k X 5 for the fixed M k and M n . In particular, A k is proportional to h 2 k \ 5 M k /M^ 
for M k M n since I(x) ~ x for x C 1. Since h rk tends to be smaller for larger values 



°This problem may also be solved by making the A5 term as an effective one through the extension of 

the model with a U(l) gauge symmetry 

f There is a severe tension between the relic abundance and the lepton flavor violation generally [TOlHTj . 

If we make neutrino Yukawa couplings small enough to suppress the lepton flavor violation, the DM relic 

abundance becomes too large. 
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Fig. 1 Contours for the branching ratio of the lepton flavor violating processes and the DM relic 
abundance in the (Mi, M v ) plane. The left and right panel corresponds to case (i) and (ii) defined 
in the text, respectively. Green dotted lines represent the contours for Br(fi —> &j) x 10 11 = 1.2,0.72 
in decreasing order. Blue dashed lines represent the ones for Br(r — > fj/y) x 10 8 = 0.68, 0.068 in 
Mi increasing order. The red solid line in the Mi < M v region corresponds to the contour for the N± 
relic abundance Cl^h 2 = 0.11 required by the WMAP. The black long dashed line represents a line for 

Mi = M^. 



of A5, A5 is expected to have values in restricted regions by taking account of the DM 
relic abundance condition as seen later. We will assume M\ < M v throughout the present 
analysis. 

The branching ratio of the lepton flavor violating process £ a — > ipj is written as [1] 

2 



3a 



6Att(G f M 2 ) 2 



k=l 



Ml 



F 2 (x) 



1 — Qx + 3x 2 + 2x 3 — 6x 2 In x 



6(1 - a;) 4 

If we use the condition (13]), we find that 



(6) 



Br(fi — > ej) 
Br{r — >■ fiy) 



3a 



6Att(G f M 2 ) 2 



0.51a 



h 2 F [ M| 

,M 2 \ , 2 ^ (M 2 



K 3 F 2 



M[ 
M 2 /A 



•(7) 



QAtt(G f M 2 ) 2 _ 

By using these formulas and eqs. (EJ), the expected values for the branching ratio of fi — > ey 
and t — > //7 can be plotted in the (Mi, M v ) plane by fixing parameters h T ±, A5, M 2 , 
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Fig. 2 Diagrams contributing to the Ni annihilation. 

and M 3 . Here we consider two cases: (i) M\ < M 2 < M 3 and (ii) Mi ~ M 2 < M 3 . In 
both cases A5 and M3 are treated as free parameters. Since h T \ and M2 are determined 
by other parameters in case (ii), this case is much constrained and predictive compared 
with case (i). 

In Fig. 1 we show the contours of these branching ratios for typical parameters. Here 
we use A5 = 6.0 x 10~ n and M 3 = 4.8 TeV. In case (i), we fix the remaining parameters 
as h T i = 1.5 and M2 = 2.8 TeV. Green dotted lines and blue dashed lines represent the 
contours of Br(fi — > ej) and Br{r — >• ^7), respectively. The former one is independent 
of Mi in both cases (i) and (ii). This is clear from the expression in eq. (J7J. Moreover, 
this branching ratio becomes sufficiently small by making M3 large enough. It should 
be noted that F 2 (Mi/M%) becomes smaller for larger M3 although larger M3 makes h r3 
larger through eq.(j5]). On the other hand, Brir — > ^7) shows different behavior in each 
case. It is independent of M v in case (ii) for the M v > Mi region. This is expected from 
the feature of which is previously remarked on eqs. ([5]) and (JTj). In case (i), Brir — > ^7) 
is not largely affected by changing M2 and M3 as long as h T i > h T 2 is satisfied, which is 
favored by the DM relic abundance as seen later. The present experimental bounds [17] 
are found to be satisfied in the wide range of parameter space shown in this figure. The 
model can be easily consistent with both the neutrino oscillation data and the bounds 
from the lepton flavor violating processes as long as parameters are suitably selected. 

Next we discuss the nature of DM in this model. Since N% is assumed to be DM, the 
condition (|3]) suggests that charged final states yielded in the DM annihilation consist 
of or r ± only. Positron and electron are only induced through the decay of these 
particles. It should be reminded again that this nature of DM is favored by the anomaly 
suggested by PAMELA and Fermi-LAT. The DM relic abundance is determined through 
the A^i annihilation, which occurs through the t-channel i] exchange diagram shown in 
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Fig. 2. The dominant contribution comes from the p-wave process. Thus, the annihilation 
cross section averaged by the spin of initial states is expressed aqf 

1 Mj^+Mj) 2 

aiv = to m^w KiAi2V ' { ) 

where we use eq. (jHJ) to derive this formula. In case (i), i\ = i 2 = 1 should be understood 
and then o\v oc h\ x . On the other hand, if the masses of N\ and N 2 are almost degenerate 
as in case (ii), coannihilation plays a role and then the contribution of = 1 and 2 should 
be summed up. As its result, we have o\v oc (h%i + h 2 2 ) 2 [13]. 

In order to estimate the freeze-out temperature Tf of N± including the coannihilation 
case, we follow the procedure given in [15], [16]. We define <r e ff and g e g as 

2 2 

i r ) 9N 1 9N 2 ^ /l±^3/2„-ij , fe 2/T /-i i r\3„ — -2x6 

Oeff = —0N 1 N 1 + 2 g < T N X N 2 I 1 + d J e + —<^N 2 N 2 (1 + e , 

5eff fleff fleff 

^cff = ^ 1 +^ 2 (l + 5) 3/2 e- :c5 , (9) 

where x = M±/T and m p i = 1.22 x 10 19 GeV. Internal degrees of freedom of iVj are 
described by g^. and 5 is defined by 5 = (M 2 — M\)/M\. If we define a e ff and 6 e fr by 
cTcfr'W = a e ff + &efr^ 2 , the thermally averaged cross section can be written as (<J e sv) = 
a e fj + 6b e s/x. In case (i), a e s and g e s are dominated by the first term since S > 0.2 [IB"] . 
On the other hand, 5 1 is assumed in case (ii). Thus, the second and third terms 
can bring the important contribution. Using these, the relic abundance of N± can be 
estimated through the formulas 

_ , 2 1.07xl0 9 x/ 0.038^ffm p iM 1 (a eff + 66 eff /x / ) 

gj m p i(GeV)(a eS + 3b eS /xf) gj xj 

where g* is the relativistic degrees of freedom at the freeze-out temperature Tf of N x . 

By using these formulas and the conditions in eq. (jSJ), we can plot the contour Q^h 2 = 
0.11 required by WMAP in the (Mi, M v ) plane. In Fig. 1, it is drawn by a red solid line 
in each case (i) and (ii) for the same values of parameters used in the estimation of 
Br(£ a — > l^i). The result in each case depends on h%x and + h 2 2 , respectively. Since 
it also depends on A 5 through the relations (jSJ) as noted before, the required fl^h 2 can be 
obtained only for rather restricted values of A5. In Fig. 1, the points on the red solid line 



g We need to remind that final states also include neutrino pairs other than the charged lepton pairs 
for the relic abundance estimation. 
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in the region satisfying both Br(fi — > e^y) < 1.2 x 10~ n and Mi < M v give the parameters 
consistent with all of the neutrino oscillation data and constraints from the lepton flavor 
violating process and the DM relic abundance. Thus, we find that the present model can 
give a very simple and consistent framework for the known experimental results. 

The values of relevant Yukawa couplings are determined for each point on the Q^i h 2 = 
0.11 line. In Fig. 1, we have, for example, 

(i) h T2 = 1.40, h T3 = 0.66 at (Mi, M v ) = (1600,1950), 

(ii) + h 2 r2 = 2.14, /i t3 = 0.66 at (Mi, M v ) = (1600, 1950). 

If we make A5 larger for the fixed Mi, larger values for M v and h T \ are required as expected 
from eq. (jHJ). On the other hand, A5 is bounded from below by the condition M v > Mi. 
Thus, neutrino Yukawa couplings are required to take rather large values by the DM relic 
abundance. This suggests that the model may be inconsistent due to these large Yukawa 
couplings, which may make the scalar potential unstable at the energy regions above a 
certain cut off scale /1. If this cut-off scale does not satisfy fi > M 3 , the present scenario 
can not work. Since larger M3 is favored from the /1 — > constraint, we can not make 
M 3 smaller enough for this instability problem. Thus, this imposes nontrivial constraint 
on the model. Since A2 is most affected by the large neutrino Yukawa couplings h T k, the 
cut-off scale \x is determined as the scale where A2 becomes negative. 

We examine this point by studying the behavior of couplings included in Cn and V to 
fix jj, using renormalization group equations (RGEs) for them. These RGEs are given in 
Appendix A. Numerical analysis is practiced for the parameters given in (ITT)) assuming 
the 0(1) values for the couplings Aj. This analysis shows that /i = O(10) TeV and then 
\i > M 3 is possible for these values of couplings at low energy regions. The scenario 
seems to be consistent with the potential instability. However, it is difficult to make \i 
much larger than M3. If N\. and rj are supposed to be suitable representations of some 
hidden non-Abelian gauge symmetry under which all the SM contents are singlet, some 
improvement may be expected for this situation. As such an example, we may consider 
SU(2) symmetry and both Nk and rj are doublet of that gauge symmetry. In such an 
extension, Cn and V are invariant and no anomaly problem occurs within these field 
contents. The RGE study of this case shows that /1 can be somewhat large. However, 
it is difficult to make \x larger than M 3 by more than one order since the running region 
of the relevant RGEs is too short. Thus, although the model can escape the instability 
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of the potential, we need to consider some fundamental model at the scale not far from 
M 3 . Since this argument on the potential stability suggests that smaller neutrino Yukawa 
couplings are favored, smaller values of Mi and A5 are also favored from the N\ relic 
abundance. On the other hand, as discussed in the next part, only a limited value of 
Mi seems to be favored from the explanation of the charged cosmic ray anomaly. This 
suggests that A5 is also required to take its value in the strictly restricted region. 

It is worthy to note that we can predict the expected values for the branching ratio 
of n — Y ej and r — > fij in this model from Fig. 1, if we can fix the value of Mi further 
by using other experimental data. Observational data on the cosmic rays from PAMELA 
and Fermi-LAT experiments may be used for such a purpose. However, if we suppose the 
PAMELA anomaly as a consequence of the annihilation of this DM, we confront difficulty 
The annihilation cross section is found to be too small to explain the PAMELA positron 
excess for the typical relative velocity of DM in the present Galaxy as mentioned before. 
In the next part we propose an extension of the model to overcome this fault. 



2.2 Extension of the model 

We consider the introduction of a new interaction which brings a large contribution to 
the Nx annihilation only at the present Galaxy and also does not modify the previously 
discussed favorable features of the model. For that purpose, we add a complex singlet 
scalar S with even parity of the Z2 symmetry!^ This singlet scalar is assumed to have 
mass Ms and couplings with other fields through the Lagrangian through the interaction 
Lagrangian 

C' N = -y k SNlP R N k - y* k S*N~ k P L N c k - M 2 S \S\ 2 - k\S\ a - (k^4> + K n rf V ) \S\ 2 . (12) 

Here we note that this is not the most general Lagrangian under the imposed symmetry. 
However, although interaction terms like <j)^<j)S and rfrjS which are not forbidden by the 
symmetry can be radiatively induced, they are largely suppressed as long as S is assumed 
to have no vacuum expectation valued In that case, S dominantly decays to N k with the 
mass M k < Ms/2. In this extended model, we find that there appears a new one- loop 
contribution to the annihilation of Ni as shown in Fig. 2. 



h The extension of the model by a singlet scalar field has been considered in other context in [15] , 
'This assumption is justified only if the tadpole diagram for S generated through the N k loop is 

cancelled by cS which can be introduced in Lagrangian. We consider such a situation here. 
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This new contribution to the Ni annihilation cross section can be estimated as 



1 M? m 2 + ml 



3 



) 



2 



{.cr 2 v) a p 



7T M% (s - M 2 ) 2 + M 2 T 2 



E 



k=l 



yiyk\h ak hf3 k M k 

{^) 2 D k 



(13) 



where the spin of initial states is averaged. We fix the final states to be charged leptons 
with masses m a ^ in this expression. This annihilation cross section is dominated by 
the contribution from the exchange of a pseudoscalar component. To obtain the total 
annihilation cross section, a and f3 should be summed up for all possible final states as 
o"2^ = YliaB^ a ^ v )aP- The definition of s, Ts and 1/D k are given by 



where r k = M%/M 2 and S is supposed to decay to the Ni pair only. This type of annihi- 
lation cross section has been suggested to be enhanced sufficiently for the explanation of 
the PAMELA data [19J. In fact, if the thermal average of (<jv) a p is estimated naively by 
replacing v 2 with a thermally averaged value - in eq. fflBl . the annihilation cross section 
shows the Breit-Wigner resonance at x r = A- through the factor [(A — ^:) 2 + 7l]~ 1 , where 
we use the definition A = 1 — -^r and 75 ~ j^lyi^A 1 / 2 . However, such a naive treatment 
has been shown to be unreliable near a resonance point [TH1 ED] ■ The enhancement of the 
annihilation cross section is overestimated in such a naive method. To obtain the correct 
enhancement, we need to calculate the thermal average 



Although this formula is derived in the center of mass system, the result is expected to 
be reliable since N\ is sufficiently non-relativistic in the present case [T6| [20] . 

In order to find the qualitative feature, it is useful to approximate this integral by 
expanding v as v = v r + v around the peak value v T = 2A 1 / 2 . Then, eqs. ( ITBl and ( fl5l) 




(14) 




(15) 
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give 



rV2 (ml + mJ)M? / * |t,,!, fc |ft„ t .ft„ t M, 



I'D 



2 

27T 1 / 2 (m 2 , + mj) /A bfc|^ a fcVM fc ^ 3/2 -xA 
- (4^ AQ /;,. U; ) X 6 ' libj 

where i/q ^ u r an d 167n/ ^ |l/i| 2 are assumed. Using this result, we roughly estimate this 
resonance effect on the annihilation cross section caused by the diagram which has N k as 
the internal fermions and r ± in final states. For that purpose, we take x r ~ 10 6 which 
is just coincident with the typical relative velocity 2 x 10~ 3 c of this DM in the present 
Galaxy. The annihilation cross section at x r is found to satisfy the relation B 

(a 2 v) ( M k \ 2 / h rk \ 4 2 



W 6 (a lV ) \M x D k J V h 



y k \\ (17) 



where h = h T \ and a/ + ti^ 2 for the annihilation and the coannihilation, respectively. 
The first two factors relevant to the masses and the neutrino Yukawa couplings of N k are 
fixed by the conditions imposed by the neutrino oscillation and the N\ relic abundance. 
The first factor is estimated to be 0(1) and decreases for larger M k . The second factor is 
considered to be less than 1 except for the coannihilation case where it can be almost 1. 
Since \yi\ is assumed small in the above discussion, we find that the desirable enhancement 
can be expected from the N 2 contribution with \y 2 \ = 0(1). These show that the sufficient 
enhancement factor to explain the PAMELA data can be obtained through the Breit- 
Wigner resonance at least in the coannihilation case. 

To obtain much quantitative estimation we calculate the thermally averaged annihila- 
tion cross section by integrating eq. f|T5l) numerically. The result for (av) = (ctiv) + (02^) 
is plotted as a function of x in Fig. 3. In this calculation we use the parameters given in 
Table 1, which can realize a point on the red line in Fig. 1. They satisfy all the neutrino 
oscillation data, the DM relic abundance required by WMAP and the constraints from 
the lepton flavor violating processes. Since the interference terms between tree diagrams 
and one-loop diagrams can be neglected in both regions Vf/c~ 0.2 and v r /c ^ 10~ 3 , the 
figure shows that we can safely use (criv) and (a 2 v) in each region, respectively. In this 



jfiere the annihilation cross section (a^) is defined as (<j 2 v) = 4(ct2w) m ± t t by taking account of all 
possible modes. See eas. (flZ2]) and (f!?51) also. 
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Fig. 3 The N\ annihilation cross section as a function of x{= 7py). The left and right panel corresponds 
to the case (i) and (ii), respectively. Parameters in the annihilation cross section are fixed to the ones 
shown in Table 1. The thin black dashed line shows the result for (av) obtained by the naive method. 

figure the result obtained by the naive method is also plotted by a thin black dashed line. 
It shows that the enhancement effect is overestimated and the annihilation cross section 
is misled to be large enough for the explanation of the PAMELA data in both cases (i) 
and (ii). However, the correct calculation shows that the enhancement can not be large 
enough for reasonable values of \yk\ in the annihilation case (i). On the other hand, in the 
coannihilation case (ii), we find that the Breit-Wigner resonance can make the annihila- 
tion cross section have a desirable value 10 -23 cm 3 /sec around v ~ v r as long as A and [ 2/fe I 
have suitable values. Here we should note that this value of A requires Ms to be finely 
tuned to Mi at the level of O(10~ 6 ). We also comment on the instability of this solution 
induced by radiative corrections. A dominant correction to the singlet scalar mass Ms at 
the one-loop level is roughly estimated as 5M| ~ j^yi^ 2 where /i is the cut-off scale of 
the model. Since [i is rather small and /i = O(10) TeV as discussed before, we find that 



Mi 


M v 


h T l 


h r2 


h T 3 


A 


\yi\ 


1 2/2 1 1 2/3 1 


(i) 1.6 


1.95 


1.5 


1.4 


0.66 


10" 6 


0.1 


2.5 0.01 


(ii) 1.6 


1.95 


0.1 


2.14 


0.66 


10~ 6 


0.015 


1.715 0.01 



Table 1 Parameter sets used to draw the annihilation cross section behavior in Fig. 2 and also to 
obtain the positron spectrum in Fig. 4. Masses are given in TeV unit. We set A5 = 6.0 x 10~ n and 
M 3 = 4.8 TeV for each case. 
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5M 2 S is the same order value as the required Ms- This means that we need the fine tuning 
of O(10 -6 ) to keep the stability of this solution from the radiative corrections. The large 
value of 1/2 may require the fine tuning up to the eight-loop order. Unfortunately, the 
model does not have any physical background to guarantee the required mass relation. 
It remains as a difficult problem how to realize this finely tuned situation from the basic 
model at high energy regions. 

It is worthwhile to stress that these values of A and \yk\ can be fixed without contra- 
dicting the required DM relic abundance which is determined by the annihilation process 
described by eq. (JE}. The reason is that different parameters are relevant to determine the 
DM relic abundance and the positron flux, respectively ^ Although the former is deter- 
mined by Mi, M v and h T k, the latter is mainly determined by and A. We note that 
the parameters relevant to the enhancement of the annihilation cross section required for 
the explanation of PAMELA and FERMI-LAT are confined to y%, y<i and Ms, although a 
lot of free parameters seems to be introduced in eq. f lT2|) . It seems to be interesting that 
these limited parameters can also allow the model to satisfy the reionization constraints 
as discussed belowjj 

It is also useful to note that it is crucial that the annihilation occurs through a one- 
loop diagram in the present enhancement mechanism. This is clear from the fact that the 
enhancement is caused by the existence of Nf., which satisfies \y).\ ^> \yi\ and h T k ~ fofc 
Although these conditions can be satisfied in the coannihilation case, in the annihilation 
case larger h T 2 requires relatively smaller h T i as seen from eq.(j5]). However, small h T \ con- 
tradicts the condition imposed by the Aq relic abundance. Thus, only the coannihilation 
case can realize fav) ~ 10~ 23 cm 3 /sec for each Mi by adjusting the values of \yu\ and 
1 2/1 1 without affecting the N\ relic abundance. 



k It is useful to note that the similar aspect is found in the case of Sommerfeld enhancements. If 

only a single annihilation channel is assumed, Sommerfeld enhancements cause a discrepancy between 

the relic density and the excesses of positron flux |21) . The present model escapes this by considering 

two processes given in Fig. 2. 

'The relevant parameters contained in eq. ([!]) have already been fixed to explain the neutrino oscillation 

data (two squared mass differences and three mixing angle), the DM relic abundance Vlh 2 and lepton 

flavor violating processes. Taking account of the supposed flavor structure, they are Mi^M^As and 

h Tl , h T2 , h Ta in the coannihilation case. 

m It is worthy to note that numerical calculation shows that {02V) has the largest value for y\ ~ 0.1 

and 0.01 at x = x r in the case (i) and (ii) respectively. It slowly decreases for larger \yi\. 
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One may worry about the potential instability caused by the large value of 1 2/2 1 shown 
in Table 1. In fact, the coupling constants k and k v can become negative at a scale smaller 
than M3 as long as 77 and Nj~ are singlets of the hidden gauge symmetry. However, if they 
are doublets of the hidden SU(2), the coupling constant y k in (a 2 v) is replaced by 2y k 
because of the gauge freedom in the one-loop diagram. This shows that a rather small 
value 2/2 — 0.86 is needed to realize the required enhancement of (cr 2 i>). We can numerically 
check that this value of 2/2 improves the above mentioned potential instability problem to 
make the extended model consistent. In this case the cut-off scale of the model is still 
determined by the behavior of A2. 

Finally we note the values of the annihilation cross section at the recombination period 

z ~ 1000, which corresponds to the DM relative velocity v/c ~ 10~ 8 . The DM annihilation 

in the period after recombination to structure formation (z ^ 6) causes the deposition of 

energy in the inter galactic medium, which brings an additional origin for the reionization 

and heating of the intergalactic gas. This additional effect is constrained from the observed 

optical depth of the universe and the measured temperature of the intergalactic gas. 

In particular, the optical depth bound brings severe constraint on the high mass DM 

as the present model since it can produce too many free electrons. If we follow the 

analysis for these constraints given in [22], the annihilation cross section should satisfy 

(02^) ~ 10 -24 cm 3 /sec for the DM with the mass 1600 GeV. In Fig. 3, we find that this 

constraint is satisfied at v/c < 10~ 5 . It corresponds to the environments in which most of 

the annihilation contribution to the relevant signal is considered to take place. Here it is 

useful to note that (c^u) does not decrease to 10~ 24 cm 3 /sec even for much smaller relative 

velocity v < v r and keep larger values than that if \y\\ ^ 0.05 is not satisfied. Thus, the 

reionization constraint rules out these cases. As long as this condition is satisfied, the 

present DM scenario can be consistent with the constraint caused by the effect on the 

reionization due to their annihilation^ In the next section we apply this extended model 

to the explanation of the anomaly suggested in PAMELA and Fermi-LAT experiments. 

"Since the Sommerfeld enhancement shows an inverse proportionality to the relative velocity of the 
two DM fields, it could cause different effects on the reionization from this model. 
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3 Positron flux and gamma ray constraints 



We estimate the positron flux yielded by the Aq annihilation following the method used 
in [23J [21] and compare it with the data obtained in the PAMELA and Fermi-LAT ex- 
periments. The positron flux in the cosmic ray at the Earth is expressed as <$> e +(E) = 
v e +f{E)/A'K (GeV • cm 2 • str • sec) -1 where v e + is positron velocity f(E) is the positron 
number density per unit energy at the Earth, which can be determined by solving the 
diffusion equation for f(E). Using the approximated solution for f(E), the positron flux 
$ e + expected from the Aq annihilation is estimated as 



*-< E » = 8^/(G*eV, E ) (t) 2 f dE ' I{ME ' E ' )] 



dE' 

where te = 10 16 sec and is the local DM density in the halo. In this study we 
use pN 1 = 0.3 GeV/cm 3 and v e + = c. Possible final states directly yielded through the 
Aq annihilation are expressed by T . dN a rp\ je + /dE' represents the spectrum of positrons 
yielded through the decay of leptons a included in the final state T . 

In this formula, the ingredients coming from astrophysics are summarized in the halo 
function I{Xd) and the positron diffusion length Xd- They are defined by 

1 (A D ) = a + ai tanh < a 2 exp + a 3 



A 2 D = AK r E | 1, (19) 

where t = log 10 (Ai)/kpc). The expressions of I(Xd) and Ad depend on the astrophysical 
model for the diffusion of positron and the halo profile [21]. In this paper we adopt med 
and isothermal profile for them to determine the parameters included in eq. f lT9|) . For 
such a model [21], parameters in A^ are K Q = 0.0112 kpc 2 /Myr and 5 = 0.70, and others 
included in I(\d) are 

a = 0.495, ai = 0.629, a 2 = 0.137, a 3 = 0.784, 

b x = 0.766, b 2 = 0.550, c x = 0.193, c 2 = 0.296. (20) 

As addressed in several work [8], the positron flux is not crucially dependent on the 
astrophysical model. We choose this model for the consistency with the constraint from 
the diffuse gamma in the cosmic ray. We will come back to this point later. 
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In eq. ( |T8|) the dependence on the assumed model for particle physics is confined in the 
factor summed up for J 7 in the E' integral. Since the annihilation cross section (o"2f) Q( g 
is proportional to m 2 + m|, the summation should be taken for 

?={e ± ,r*), (ftV), (t+,t-), (21) 

which can yield positrons finally. This feature is caused by the flavor structure of neutrino 
Yukawa couplings Since smaller h r3 and larger M 3 are favored from the fi — > e^y 
constraint, the N 3 contribution to the loop effect may be neglected. If we take account of 
these and also assume that \y 3 \ is sufficiently smalj°|, the positron flux $ e + due to the Ni 
annihilation can be expressed as 

$ e+ ~ 1.25 x 10- 3 (a 2 v) 



AO 2 GeVV (\ TeVY 

v e ) v^T; 



Mi 

dE' I(E, E') 



ldiV |t + je + + 3 fdN T +^ e + dN T - 



(22) 



4 d£' 4 V dE' dE' 

where (GeV • cm 2 ■ str • sec) -1 is used for the unit of $ e + and the total cross section 
((j2v){= 4(cx At ± T T)) is determined in our extended model as 

16 (GeV) 2 Mfml ff M^ V 



( 47r ) 5 (A-^) 2 +7 2M , 4M Hfe D * 



Here it should be noted that we can keep the favorable feature such that the final states 
of the Ni annihilation consist of heavier leptons onlyj^l The fact that e ± are not directly 
produced is favored to explain the Fermi-LAT data, which show no bump in the hard 
e + + e~ spectrum. The directly produced tend to be much harder than indirectly 
produced energetically. The positron spectrum given by (122]) has large contributions 
from r + decay, which causes a softer spectrum for the final positron and electron spectrum. 
The decay of ^ to yields much harder positron than the decay. The concrete 
model with these mixed final states seems not to have been considered in the analysis of 
the anomaly suggested through the PAMELA and Fermi-LAT experiments, 



11^ 





°Under this assumption, we can safely neglect the ./V3 contribution to the one-loop annihilation dia- 
gram. In this case N\ annihilation does not yield positrons directly. 

p Although higher order radiative corrections can induce coupling of S with the ordinary Higgs scalars, 

their effect is small enough to neglect them in the analysis of Ni annihilation. 

q Although model independent analysis for this kind of mixed final states is found in the paper by 

Meade et al. [5] (see Fig. 12 in it), it is not based on a concrete particle physics model. 
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Fig. 4 Left and right panels show the predicted positron excess at the PAMELA regions and the 
prediction for the e + + e~ flux at the observation regions of Fermi-LAT and H.E.S.S., respectively. In 
both panels, DM mass Mi and annihilation cross section (atv) are fixed as (Mi (TeV), (<T2v) (cm 3 /sec)) = 
(1.2, 4.1 x 1(T 23 ), (1.6, 6.7 x 1CT 23 ), (2.0, 9.2 x 1(T 23 ). The normalization of background fluxes is taken 
to be N<j, = 0.64. 

The energy spectrum of positron — r£r~ can be computed by using the PYTHIA 

dE 

Monte Carlo code [25] . We determine the positron spectrum by fitting these simulation 
data for both the fi + and r ± cases. Details of the analysis are given in Appendix B. We 
apply this result to eq. ( 1221 to find the positron flux $ e +. We fix parameters included 
in the cross section (<72f ) by using the ones which realize the point in the allowed region 
shown in Fig. 1. They are also summarized in Table 1 for the case of M\ = 1.6 TeV. 
As expected background fluxes for positrons and electrons, we use the empirical formulas 
given in [26J, 

4.5£ a7 



ibkg _ "' - LU - L7 I T\T "' ' 

e~ , t -i pn.<) i QOP2.15 T $ 1 i nnci1.fi i cnn i conc'4.2' ^ > 



i + 650£ 2 - 3 + i500£ 4 - 2 ' 

0.16E- 11 0.70E - 7 
1 + HE - 9 + 3.2E 2 - 15 + *1 + HOE 1 - 5 + 600E 2 - 9 + 580E 4 - : 



where E should be understood in a unit of GeV and iV$ is a normalization factor. 

Using these formulas, we plot the positron fraction $ e +/($ e + + $ e -) and the total 
flux of e + + e~ scaled by E 3 in two panels of Fig. 4, respectively. In the left panel, the 
data points for positron excess of PAMELA [2], CAPRICE94 [28] and HEAT95 [29] are 
also plotted. On the other hand, the data points for the e + + e~ flux of Fermi-LAT [3] 
and HESS [27] are plotted in the right panel. In this figure {(T^v) is fixed to make the 
positron flux $ e + to realize a good fit to the data of PAMELA and Fermi-LAT for each 
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Mi value. This figure shows that the flux of positrons and electrons predicted from the 
annihilation of iVj. in this extended model can give rather good fits with these experimental 
data. Especially, the predicted flux fits well both data of the PAMELA and Fermi-LAT 
experiments for Mi = 1.6 TeV and (a^u) = 6.7 x 10 -23 cm 3 /sec. This value of (cr^v) can 
be realized by the parameter sets in case (ii) given in Table 1. If we apply this information 
for Mi to Fig. 1, we can predict the value of Br(fi — > cy) and Br(r — > fi^f). The figure 
shows that the predicted value for fi — > ey is within the reach of the MEG experiment 
[30] . Thus, lepton flavor violating processes could be a crucial probe for this model. 

Annihilation of the DM can cause additional contributions to the cosmic gamma ray. 
In fact, if hard charged leptons are produced as the final states of the DM annihilation, 
high energy photons are also produced through several processes. One of their origin is 
the inverse Compton scattering of positrons with CMB, star light and interstellar photon 
[3T] . The other ones are final state radiation or internal radiation [33]. The gamma ray 
flux expected from the former one does not depend on the particle physics model as long 
as the positron flux data presented by PAMELA is assumed. It can be used as a crucial 
constraint on the model. Since the gamma ray flux caused by the latter ones depends 
on the adopted particle physics model, the predicted photon spectrum can be used to 
discriminate the model from others on the basis of the deviation of the photon spectrum 
from the one of background in the future observation. 

In the present scenario, the DM has mass of 0(1) TeV and it can decay into t^. Thus, 
substantial constraints are expected to be imposed by the gamma from the former origin 
and also the gamma produced through the decay of tt° which comes from the decay. 
These give strong constraints on the gamma ray flux at higher energy regions. Various 
studies related to this issue have been done in the model independent way or in the fixed 
models [3TJ |3U [22] . The constraints obtained from analyses of the first year of Fermi 7-ray 
observations are also given in [35]. Their results for the gamma ray flux associated with 
the DM annihilation into charged lepton pairs are applicable to our model to examine the 
consistency with the diffuse gamma ray observations. They show that the galactic diffuse 
gamma data constrain the assumed DM halo profile severely. Only the restricted halo 
profile called isothermal seems to be consistent with the observations. In fact, following 
the study by Cirelli et al. in [35] for the cases with the final states and r + r~ for the 

DM with mass around 1.6 TeV, the DM annihilation cross section (cry ) is shown to be less 
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than 6 x 1CT 23 cm 3 /sec and 1 x 10~ 22 cm 3 /sec, respectively. Papucci et al. in [35] gives 
much stronger constraints for the r + r~ case. The results in this section show that our 
model can satisfy these constraints for / u + but the situation seems to be marginal for 
the r + r~ case. Thus, the present model may be considered to work well in the isothermal 
profile, although this type of halo profile is considered to be disfavored by the iV-body 
simulation. We also note that the diffuse neutrino flux satisfies the present observational 
constraints [311 E2] ■ 

4 Summary 

The radiative seesaw model is a simple and interesting extension of the SM by an inert 
doublet scalar and singlet fermions. It can give the origin of both small neutrino masses 
and DM consistently. However, if we try to explain the positron excess observed by the 
PAMELA experiment on the basis of the DM annihilation in this model, an extremely 
large boost factor for the annihilation cross section is required. In this paper we have 
proposed a simple extension of the model by introducing a singlet scalar. In this extended 
model, the DM annihilation cross section can be enhanced in the present Galaxy through 
the Breit-Wigner resonance without disturbing the features in the original model, which 
are favored by the neutrino masses, the lepton flavor violating processes and the DM 
relic abundance. However, it should be noted that the mass of the singlet scalar has 
to be finely tuned at the level of O(10 -6 ) for this enhancement. Final states of DM 
annihilation are composed of heavier leptons only and the ratio of fi + and r^ 1 contribution 
to the annihilation cross section is 1 to 3. As a result of these features, the data for the 
positron and electron flux observed by PAMELA and Fermi-LAT are well explained in 
this extended model as long as the coannihilation among Nf. occurs. It is interesting 
that these results are closely related to the flavor structure of neutrino Yukawa couplings, 
which induces tri-bimaximal mixing. 

This extended model may be checked through the study of lepton flavor violating 
processes such as fi — > ej in the MEG experiment and others in near future. The cosmic 
positron and electron flux at higher energy regions may be clarified by the future CALET 
experiment, which can observe e ± flux up to 10 TeV [36]. Viability of the model may 
also be confirmed through this experiment. Although the diffuse gamma ray flux imposes 



20 



severe constraints on the model, they could be consistent as long as the specific halo 
density profile called the isothermal profile is assumed. Detailed knowledge on the density 
profile of the DM halo seems to be required to judge the validity of the explanation given 
here for the anomaly reported by PAMELA and Fermi-LAT. 

This work is partially supported by a Grant-in-Aid for Scientific Research (C) from 
Japan Society for Promotion of Science (No. 21540262). Numerical computation was par- 
tially carried out by using the computing facility at Yukawa Institute. 
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Appendix A 



In order to study the stability of the scalar potential, we need the renormalization 
group equations (RGEs) for the coupling constants included in the scalar potential V. 
This model is an extension of the SM with three singlet fermions N*. and one doublet 
scalar 77. Thus, the RGEs are similar to the ones of the ordinary two doublet Higgs 
model. However, we need to take account of the effect of large neutrino Yukawa couplings 
h a k, which are assumed to satisfy the relation (j3j). 

If we assume the existence of the hidden sector gauge interaction mentioned in the text, 
invariance of Cn restricts the representation of Nk and i] to be an adjoint representation 
of SU(iV) or a doublet of SU(2), for example. In case of the adjoint representation, scalar 
potential V should be modified. No anomaly problem appears in both cases. However, 
if we note that one-loop diagrams with internal lines of and r\ have additional group 
theoretical factor dim(i?), we find that the latter is favored from the constraints of lepton 
flavor violating processes. 

To prepare RGEs applicable to this extended situation, we assume that Nk and r\ are 
singlets (N = 1) or doublets (N = 2) of a hidden gauge symmetry SU(2). All SM fields 
are singlets under this group. A set of relevant RGEs can be written in the following form 

16tt 2 
16tt 2 



~dT 



24A 2 + 2NX\ + NX\ + 2iVA 3 A 4 + 12Ai/i 2 - Qhf + k 2 



8(N + 2)A 2 + 2A 2 + A 2 + 2A 3 A 4 + 4A 2 



2[Ki + K2 + 



3C 2 (N)gl 



2N-2) 



-9 h + < 



16vr 2 ^ = 4(3AiA 3 + (2A^ + l)A2A3 + AiA4 + A^A 2 A 4 ) + 4A 2 . + 2A 



16tH 



167T 



16tt z 



, <iA 4 
~dt 

,dn 
~dt 

dn^ 
dt 



+2A 3 
2A 4 



3h 2 + 2 f + h 2 r2 + h 2 r3 



3C 2 (N)g 2 h 



2 (Ai + A 2 + 2A 3 + A 4 ) + 3h 2 t + 2 [ h 2 Tl + h 2 r2 + -h 2 3 ) - 3C 2 (N)g 2 h 



+ 2k^,k^ 1 
3 



20k 2 + 2k% + 2Nk 2 v + 8N ^(^M 2 - 4|y fc | 4 ), 



k=l 



+ ±Nk v X 3 + 2Nk v X 4 + Skk^ + 4kJ + ANk^ ^ \y k \ 2 + 6/^/1 



k=l 
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16tt 2 ^ = (87V + 4)k„A 2 + 4re A 3 + 4k A 4 + Skk v + Ak 2 + Ak v (h 2 rl + h 2 r2 + ^/i 2 3 ) 



+47V S l^l 2 - 32 ( h lM 2 + KM 2 + 2^ 3 ll/3| 2 ) - QC 2 (N)k v9 1 

k=l 

16tt 2 ^ = / i „((7V + 4)(^ 1 + ^ 2 ) + 3^ 3 + 2|^| 2 -3C 2 (7V)^) (i = l,2), 
16vr 2 ^f = ^3(2^ + ^ + ^(7V + 4)^ + 2|y 3 | 2 -3C 2 (7V)^ 

16tt 2 ^ = y i (8\y l \ 2 + 2NY,\yk\ 2 + 2h 2 Tl (2 + 5 3i )^ (i = 1,2,3), 



IGtt 2 ^ = ^(-117V+^2T(7V)) ; 



JVfe,»7 



16vr 2 ^ = -7, 3 3 , (25) 

where 62(71) and T(n) stand for values of the second order Casimir operators defined by 
Y ja T a T a = C 2 (n)l and tr(T a T b ) = T(n)5 ab for SU(n) generators T a in the fundamental 
representation. Thus, = for the TV = 1 case, and C 2 (2) = | and T(2) = | for 
the N = 2 case, respectively. In these RGEs we take account of the contributions to 
(3- functions only from the top Yukawa coupling h t , the neutrino Yukawa coupling h a k, the 
strong gauge coupling g 3 and the hidden gauge coupling except for the couplings in 
the scalar potential. Since the /3-function of A5 is proportional to A5 due to the symmetry 
discussed in the text, it is kept sufficiently small to be neglected. 

Appendix B 

In the present model the final state positron is yielded as a consequence of /i + and 
decay. We determine the energy spectrum of such positrons by using the PYTHIA Monte 
Carlo code [25] . If we write an expectation value of the number of this yielded positron 
per the decay of a(= y u + ,r ± ) as N ae +, PYTHIA gives the positron spectrum 



dE 

The spectrum obtained from this simulation is shown in Fig. 5, where the result for the 
a = fi + is plotted in the left panel and the one for a = pair is plotted in the right 
panel. We find from these figures that the positron produced through the decay of r 1 * 1 is 
softer than the one for fi + as mentioned in the text. 

In order to fix their empirical formulas approximately, each data set in Fig. 5 are fitted 
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Fig. 5 The energy spectrum obtained for the DM mass M\ — 1.2, 1.6, 2.0 TeV by simulation. 

dE 

Left and right panels show the positron spectrum obtained from the decay of the n + and pair, 
respectively 



by using the functions 



dN a „+ 



E 



d n (M 1 - E) 1 / 2 



(26) 



dE ^ (E + E r 

where E is a constant and E should be understood in a GeV unit. As results of this 
fitting, we find that the coefficients d n in the above fitting functions should take the 
values shown in Table 2. We have A^+ )6 + = 1, N T ± <e + ~ 1.3 by integrating the obtained 
spectrum. This corresponds to the fact that the decay of r + is composed of various modes 
such as r + — > e + u T u e , — > hadrons — > e e e T , while the decay mode is dominated only 
by /i + e+v^Ve for 



particle(ct) 


M+ 


r± 


Mi (TeV) 


1.2 


1.6 


2.0 


1.2 


1.6 


2.0 


do 


-2.76 x 10" 3 


-3.86 x 10" 3 


-2.13 x 10" 3 


-2.50 x 10" 6 


-1.14 x 10" 6 


-7.75 x 10~ 7 


di 


4.13 x 10 1 


1.20 x 10 2 


7.19 x 10 1 


2.94 x 10~ 3 


1.85 x 10" 3 


1.44 x 10" 3 


d 2 


-1.52 x 10 5 


-9.26 x 10 5 


-6.01 x 10 5 


2.43 x 10° 


3.27 x 10° 


3.86 x 10° 


E 


7.12 x 10 3 


1.51 x 10 4 


1.64 x 10 4 


6.54 x 10 1 


9.60 x 10 1 


1.24 x 10 2 



Table. 2 The coefficients d n determined by the fitting to Fig. 5. 



In the text we make the estimation of the positron flux (1221) by using these positron 
spectra. The parameters included in the cross section (o"2f ) is determined by a point in 
the allowed regions (on the red solid line) shown in Fig. 1. Other parameters A and yu 
relevant only to the N\ annihilation at the present Galaxy are fixed to make (p-iv ) a suit- 
able value O(10 -23 ) cm 3 /sec for the explanation of PAMELA data. They are summarized 
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Table 1 in the case of Mi = 1.6 TeV. 
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